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SUMMARY 


The governing differential equations of motion for a high speed 
cylindrical roller bearing are developed under the assumptions that 
the bearing is isothermal and that the roller tilt and skew are very 
small. Two sets of differential equations are presented; the first 
which deals with planar motion of the roller bearing system, and the 
second which includes the effect of roller skewing. The equations 
as presented are in a format for programming on a digital computer. 



INTRODUCTION 


Current design trends for aircraft gas turbine engines are to 
achieve higher speeds and lighter weights. The primary advantages 
of these higher speeds and light weights are to increase payload 
and to reduce the drag induced by the externally mounted machines. 

The mainshaft bearings used for high speed turbines are both 
cylindrical roller and ball bearings. The roller bearings have 
both high radial load capacity and the ability to permit rela- 
tively unconstrained axial motion under differential thermal expan- 
sion of the rotor when the engine is at operating condition. At 
high speeds and under low loads the rollers in the bearing have a 
tendency to slip or skew or both. These phenomena may cause the 
bearing to fail prematurely due to surface distress and/or jamming 
between the cage and the roller which can lead to catastrophic 
bearing failure. 

In lightly loaded roller bearings operating at very high 
speeds, the rollers are subjected to centrifugal forces which keep 
them in contact with the outer race. The rollers are also sub- 
jected to contact forces as they traverse through the load zone of 
the bearing and react against the externally applied loads. These 
loading conditions, together with the effects of elastohydrod3mamic 
lubrication, influence the rollers to slip relative to the races. 

As a consequence, the cage rotational speed falls below the theo- 
retical epicyclic speed and, as the rollers traverse through the 
load zone, they are subjected to a varying radial force which 
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causes them to speed up and to slow down. 

To understand the behavior of high speed lightly loaded roller 
bearings, it is necessary to examine the dynamic behavior of the 
rollers and the cage, and to determine the force interactions 
present in the bearing. Methods which are currently available are 
steady state force and moment balances. These analyses predict 
the average forces and moments on the rollers and cage. Transient 
effects such as cage whirl, roller slip as a function of angular 
position, or skewing cannot be predicted with any reasonable accu- 
racy using these steady state methods. A time-domain solution to 
the differential equations of motion of the rollers and the cage 
in the bearing will provide detailed information on the cage-roller 
and roller- race interactions in the bearing. This detailed knowl- 
edge of the behavior of high speed roller bearings will aid engi- 
neers and engine designers in finding ways to improve their design 
with the objective of increasing roller bearing operating life and 
reliability. 

Two articles have reported studies of the transient 
dynamic behavior of ball bearings. Walters [1] investigated the 
interaction between the cage and the ball while Gupta [2] investi- 
gated the ball- race dynamics without the effects of the bearing 
cage. Rurabarger, et al. [3] modeled the roller bearing and per- 
formed a steady state dynamic analysis on the resulting system of 
equations. The results of this study showed that the predicted 
cage slip compared favorably with experimental data on cage slip 
when the effective viscous drag forces of the lubricant on the 
cage were adjusted by varying the effective density of the lubricant. 
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Recently, Gupta [4] published a series of papers on the dynamic 
analysis of rolling element bearings. The approach to the problem 
utilized by Gupta is essentially the same as the approach utilized 
in this report. The work by both investigators was conducted inde- 
pendently and resulted in the same general approach to the problem. 
Most of the models in this report and in [4] are based on the work 
of Rumbarger et. al, [3], 

The major assumption used in this work is that the bearing 
system is at a thermal steady state. This implies that the lubri- 
cant properties do not vary because of temperature changes and that 
the geometry of the rollers and races remains constant because of 
no thermal expansion. The motion of the rollers was assumed to be 
essentially planar. The tilt and skew angles of the roller are so 
small that they do not appreciably change the dynamics of motion 
of the roller, but the tilt and skew angles do have a dominant effect 
on the forces acting on the roller. 

The purpose of this program of research was to achieve the 
following goals: 

(1) To develop models which describe the interactions among 
forces, relative displacements and relative velocities 
between components in a cylindrical roller bearing and 
which are computationally efficient, 

(2) To develop a working computer program for planar motion 
of the rollers and cage in the bearing which incorporates 
the models from goal (1), 

(3) To develop the additional models and governing differen- 
tial equations which describe the three dimensional 
motions of the rollers. 
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(4) To develop a working computer program for the three dimen- 
sional motion of the components of the bearing. 

The above goals were met and the results are presented in the 
subsequent sections of this report. The governing relationships 
between forces and displacement and velocity for the many conditions 
present in roller bearings are presented in the section THE CONSTI- 
TUENT MODELS. 

The governing differential equations of motion are developed 
for the two dimensional motion. Then the governing differential 
equations of motion for the three dimensional motion are derived. 

The assumptions upon which these equations are based together with 
the process of non dimensionalizing the equations are discussed. 

The final section of this Part contains the conclusions and 
recommendations for further work. Part II [5] of this report con- 
tains the organization of the computer program and results obtained 
by using this program. 
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THE CONSTITUENT MODELS 


The primary forces present in a high speed roller bearing are 
the normal forces between the roller and the races needed to balance 
the external forces and moments and the centrifugal forces on the 
rollers. Associated with these normal forces are the tractive 
forces between the rollers and the races. 

The cage-roller interaction forces result from a lubricant 
squeeze film phenomenon and, in extreme cases, a contact phenomenon. 
The remaining forces are a result of shear stresses on the surfaces 
of the cage and rollers due to a shearing of the fluid medium. 

The relationships between the forces and the deflections and 
relative velocities are based on the global coordinate system shown 
in Figure 1. The y axis is placed in a vertical position, parallel 
to the gravity vector, and z axis is parallel to the centerline of 
the shaft. The x axis is perpendicular to both the y and z axes. 

The angular velocities of the j ' th roller and the cage are shown in 
Figure 2; ^ ^ is the angular velocity of the roller about its own 

axis, ^ . is the angular velocity of the center of mass of the roller 
about the global z axis and is the angular velocity of the cage 
about the global z axis. The forces and moments acting on a roller 
are shown in Figures 3 and 4. Figure 3 shows the net forces and 
torques projected on a plane normal to the z axis while Figure 4 
shows the distribution of load between the roller and inner race 
and the roller and outer race together with the roller-guide flange 
impact forces (top view only) . The net forces and torques acting on 
the cage are shown schematically in Figure 5. The net translational 
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Roller-Guide Flange 
Impact Forces (4) 


Traction Force Distribution 
on the Outer Race 


Load Distribution on 
Outer Race 


Slices of the Roller 


Load Distribution on 
Inner Race 


Figure 4. Roller-Race and Roller - Guide Flange 
Interaction Forces 





forces, and ^Fy, are the sum of the forces acting on the cage 

as a result of cage-roller interaction and cage-guide flange inter- 
action, The torques and T__. _ are the sums of all the driving 

DKlVlif DRAG 

torques and drag torques respectively as a result of cage-roller inter- 
action, cage-guide flange interaction and drag torques induced by the 
cage moving through the lubricant/air medium. 

This section identifies each of these forces and presents the 
detailed relationships necessary to calculate them. 


Load distribution in the rollers 

When the spacing between the inner race and the outer race is 
less than the undeformed diameter of the roller, a normal contact 
pressure between the roller and the races is induced. The movement 
of the inner ring relative to the outer ring may be a result of 
directly applied external forces or overturning or misalignment 
moments. The contact pressure on the rollers will be sjrmmetrical 
if the loading is due only to a directly applied radial force. If 
misalignment is present, the contact pressures will be non- symmetrical 
along the length of the roller. In this investigation only crowned 
rollers are considered. 

This problem was investigated by Harris [6] and by Liu [7] who 
considered the effect of centrifugal force for high speed cylindrical 
roller bearings. The roller is divided into segments or "slices" 
along the axis and the contact pressures are determined over each of 
these segments. Figure 4. This provides an approximation to the non- 
uniform load distribution along the axis of a roller under misalign- 
ment. 
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The analysis is based on the work of Palmgren [8] which gives 
a relationship between the deformation (g) of a cylinder of length 
(jj) and an applied load (Q) as 


g = A 




0.9 


0.8 


( 1 ) 


where A = 5.46 X 10 for steel components. As the computer program 
will utilize displacements as the state variables, the above expres- 
sion is inverted, utilizing a load intensity q. 



( 2 ) 


then 


1.11 

g 

1.11 0.11 

A i 


(3) 


When the bearing is subjected to radial displacements and 
^ and angular misalignments 0^ and 0^, the equations (13) and 
(14) in Liu [7] are replaced with the following expressions for 


the total contact deformations at the inner and outer contacts 
respectively at the center of the i*"^ slice of the roller. 
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Ay cos - ■^/2) 

+ fz- + R (1/2 (0 cos d> . - Q sin^.))l « 

L 1 y X j J 

cos<f>^ - 0^ “ 3j (cos (<^J +^/2 - §^))J 


- - c. 


J 1 


( 4 ) 


S2ij “ +R <1/2 (0y cos<^^ - 0^ sinV^^))] 


Pj sign (cos + ^/2 ' §o))] 


- C 


(5) 


where is the initial separation due to crowning 
and is the diametral clearance 

and where is given by the following relationship based on Ref . 5; 



(6) 


if 0 =0 and 0 > 0 , $ = ^/2 

0=0 and 0 < 0 $ =-'^/2 

X y » lo " 


Thus, using the conventions in Liu [7], the misalignment 
applied to any roller at any angular location is 


0. = © cos <3^, + 0 cos (d^, + ^/2) 
J y 'J X ^rj ' ' 


(7) 
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The centrifugal force on a roller Is given as (based on equation 
(17) in Liu [ 7 ]) 


NP 


= "/8 ("/30)2p <2(R + r))o)^^ (2r - 2C^)' 


k=l 


( 8 ) 


where p is the mass density of the roller 

is the angular velocity of the cage 

is' the roller crown drop at segment k as defined for 

equation (4) 

The load distribution on the rollers in the bearing is deter- 
mined by enforcing the conditions of static equilibrium of forces 
and moments on each of the rollers. The diametral clearance, 
initial separations due to crowning and linear and angular defor- 


mations 


of the inner race with respect to the outer race (A^> Z\yi 


0,0) are given as inputs to the calculation. The computational 
y o 

procedure will give the central deformation of the roller, 

and the angular displacement of the roller, for all the rollers. 

The contact pressure distributions are then calculated using the 

and the load and moment on the bearing can be 

determined from equilibrium considerations of the inner contact 

forces on the inner race. 

The static equilibrium of force on the roller implies that 
the sum of the forces at the inner contact and the centrifugal 
force must balance the forces at the outer contact. Setting these 
sums equal to a residual (^j) gives 
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= z r — 

J i=l ,1.11 0. 
u A 0 


11 


(glj_j) 


1.11 


, 1.11 0.11 

A I 


f , 1.11 

J 


(9) 


+ F 


j =1, ..., NZ 


The static equilibrium of moments on the roller implies that 
the sum of the moments of the inner contact forces must balance the 
sum of the moments of the outer contact forces. Setting these sums 


equal to a residual gives 


NP 


F 

j+NZ 


Y w z. , 

• 1 . 1.11 . 0.11 
1=1 L A i 


11 


W iCf ^ 

1.11 0.11 


1.11 


( 10 ) 


j = 1, ..., NZ 

where NZ is the number of rollers 

NP is the number of slices along the axis of the roller 

w is the width of a slice = — 

NP 

When the state of static equilibrium is reached, the residuals 
will be reduced to zero. This is a system of (2 X NZ) nonlinear 
algebraic equations in the NZ variables and the NZ variables 

Defining 


r w 

* 1.11 .0. 


11 


( 11 ) 
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equations <(9) and (10) become 


NP 

riVr, .1.11 . xi*ii1 


1 . 11 - 


+ F 


i=l 


j ” 1» NZ 


NP 




i=l 


( 12 ) 


(13) 


j 


1 , 


NZ 


A Nevjton-Raphson procedure is used to solve this system of 
equations. The derivatives of the residual functions to be used 
in the Jacobian matrix are given as: 


SF 

-J— =0 if U4=j 




(14) 


NP 


i = l 

. I r z • "*■ “I" (© <j> . - Q sin cf> .)~\ 

iLi 2^y r 2 X 'J 

sign (cos (^^ +^l2 - 


if k=j 

4— = 0 if k=^j 




( 15 ) 


NP 

i=l 


if k=j 
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0 


(16) 


S3k 


= 1 


if k4j 

NP 

.11 r - z, (g,ij)"-“} 

'{ [^1 "T % ■ ®x 


sign (cos(^ ■'■ " 2 " ■ ^ ))} 1^ 


B F . + NZ 

J = 0 if k4j 

3 5; 


( 17 ) 


NP 

1 11 pY r , ^0.11 . ^o.ll1 , . 

= 1.11/ L.-z, (St,,) - -* 


i=l 


'i ^^lij' 


Denoting changes in the values of and by S 3j 

respectively, the system of equations to be solved is of the form: 



r 1 

r 1 


— -1 


1 1 
1 1 

I_J 



_ 3 9 k J 







^''j+NZ 




_ Ti. 

1 


— I 

CD. 

/D 

L 3 55 kl 



1_- j-wzjj 


(18) 


The submatrices in the coefficient matrix are diagonal, so 
the solution of the above system of 2 X NZ equations in 2 X NZ 
unknowns can be reduced to the solution of NZ (2 X 2) systems of 
equations, A solution may be obtained very efficiently. This 
procedure is programmed in subroutine ROLDIS [5] and its sub- 
sidiary subroutines. 
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The total normal forces on the roller at the inner race and 


the outer race are then given as: 


N. 

in 


NP 

(j) - V _w 


■ i-fc 


^ 11 0.11 ^%ij^ 

i=l ^ ^ 


1.11 


NP 


N 


out 


y , 1.11 . 0.11 


A A 

x-1 


( 19 ) 


j = 1 , . . , , NZ 

T raction between the rollers and the races 

The tractive forces in an EHD contact between a roller and a 
race are dependent on the lubricant formulation, the Hertzian contact 
pressure, the rolling velocity, the sliding velocity and the temper- 
ature of the surfaces. Experimental evidence, e.g. [9], has shown 
that the tractive forces are not directly proportional to the shear 
rate of the lubricant over a wide range of shearing rates. For a 
constant maximum Hertzian pressure and rolling velocity, the traction 
force rises, peaks and then decreases with increasing sliding velocity 
An approximate model which was utilized for the development of 
the roller bearing dynamics computer program is based on a curve of 
normalized traction coefficient versus normalized sliding speed. 

This curve was presented in Ref. [10] for a Mil L-7808 oil. 

The model is a piecewise linear approximation to the above ref- 
erenced curve which is presented in Fig. 6. The normalizing constants 
and u^ are given as [10]: 
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Relative Traction Coefficient /i/// 



0 12 3 

Relative Sliding Rate 

Us/Us 


Figure 6. Normalized Traction Curve for Mil-L-7808 Oil 
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( 20 ) 



4.598123 

fCq^) = 0.802 (q^ X 10"“) 


where q 


u 

s 


* 

u 

s 

* 

u 

"b 


V 


Vo 


* 

M- 


* 


2 

is the maximum Hertz contact pressure (n/m ) 

is the sliding velocity at a roller/race contact 
(m/s) 

is the normalizing value of u^ 

is the base value of the relative sliding velocity 

upon which the correlation is based. (For Mil 

L-7808 oil, u = 1.397 m/s 
®b 

is the rolling velocity at the roller/race contact 
(ra/ s) 

is the oil viscosity at operating temperature 
(m^s) 

is the traction coefficient 

is the maximum value of the traction coefficient 
corresponding to u^ 

is a base value of the traction coefficient upon 
which the correlation is based. (For Mil L-7808 
oil, = 0.019). 


( 21 ) 


( 22 ) 
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The traction coefficient is given by the following model: 


0 



^ 0.25 



2.80 


u 

s 


u * 
s 


u 

0.25 < —V ^ 0-50 

'V. U * ^ 


— = 0.70 + 0.80 

\i* 



0.25 


) 


u 

0.50 < — V ^ 0.75 
s 


— = 0.90 + 0.30 
M.* 




u 

0.75 < — ^ < 1.00 
s 


II 


= 0.975 + 0.10 



0.75 


) 


1.00 < 


u 

s 


u * 
s 


"il* 


1.00 


0.0571 



1.00 


) 


(23) 


u 

s 


u * 
s 


> 15 


pi* 


1.00 - 0.0571 ^15.0 - l.o'^ 
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Although this model is an approximation, it fairly represents 
the traction phenomenon in an EHD contact V7ith Mil L-7808 oil. 

Other models can be utilized in the roller bearing dynamics program 
for different oils merely by changing a subroutine. The above equa- 
tions are programmed in the function subroutine COEFF [5]. 

The traction force between the rollers and the race is of the 

form 


F 

IN 


U.. N. 
f^in in 


F 

OUT 


M-out 


out 


(24) 


Drag forces on the cage 

As the cage rotates in the bearing it is subjected to forces 
which accelerate its motion and to forces which retard its motion. 
This section is devoted to the description of those forces xi?hich 
retard the motion of the cage excluding those caused by cage/roller 
interaction. . . * 

These drag forces are induced by the shearing of the air/oil 
mist mixture enveloping the cage of a high speed bearing. They 
are considered to act on the outer and the side surfaces of the 
cage. 

Boness [11] experimentally described and analytically discussed 
the influence of the amount of lubricant supply on cage and roller 
motion in high speed roller bearings. Poplawski [12] introduced a 
fluid churning loss and cage pilot surface friction formulations. 
Rumbarger et. al. [3] extended Poplawski 's approach and further 
examined the influence of fluid drag forces on individual surfaces 
of the cage. 
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The models used in this computer program are taken from Rum- 
bar ger et, al. [3]. The drag torques acting on the outer cylin- 
drical surface of the cage are described by an equation of the form 


^CDO *^7 ^ ^OUT 


(25) 


where is the drag torque on the surface 


w 


OUT 


is the wall shear stress 

is the surface area 

is the outer radius of the cage 


The wall stress for a surface rotating in a viscous fluid is 
defined by [13], 

2 


_ f p u 


w 


( 26 ) 


where f is a friction factor 

/ 

p is the fluid density 

u is the mass average velocity of the fluid 


and u = 


^OJT^ 


for this correlation. 


where 'd' is the angular velocity of the cage. 

The effective density of the oil/air mixture about the cage 

is given in [3] as 

2 


-= P (DECFUL) 

*^eff (0,4 + 0,6 DECFUL) 


(27) 


where DECFUL is the ratio of the oil volume in the bearing to the 
total volume of the bearing, 
p is the density of the oil 

is the effective density 
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The Couette turbulence regime is assumed for the outer surface, thus 


where 




0.85596 


N, 


RE 


r ujc 
c 

V 


and 



Here 



is the Reynolds number 

is the laminar friction factor 


(28) 


(29) 

(30) 


If the Reynolds number is less than 2500 or the Taylor number 
is less than 41, the appropriate friction factor is taken to be the 
laminar friction factor f^. The Taylor number is defined as 


r a)C 

N = —2^ 

TA u 


VT 


(31) 


where r is a characteristic radius , equal to r 

c 

(U is the angular velocity, equal to 

c is a radial clearance 

V is the oil kinematic viscosity 


OUT 


For the side walls of cages, the characteristic radius to be used 
in the correlations below is modified by using the relationship 


5 _ 3 . 2 2x 

’^c ^OUT ^^OUT " ^IN ^ 


( 32 ) 


where ^quT outer radius 


r,,, is the inner, radius 
IN 
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The Reynolds number is given in this case as; 


RE 


( 33 ) 


and the drag torque on the sides is given as 


where 


^CDS = -T S 


(34) 



3.87 




T7T 


0.146 


<^RE> 


T7T 


for laminar flow, 
Nre < 300,000 


for turbulent flow 
\e ^ 300,000 


(35) 


These relationships are programmed in subroutine CGTRQ [5]. 
Roller- cage interaction forces 

In roller bearings the rollers have the freedom to move relative 
to the cage. In the load zone, the roller may be accelerated and 
approach the front part of the pocket resulting in a driving force 
on the cage. Outside of the load zone the roller may start to slip, 
decelerate and then approach the rear part of the pocket resulting 
in a retarding force on the cage. In many instances the squeeze 
film forces in the lubricant which are generated by the relative 
approach of the roller and the cage are not large because of the 
relative freedom of motion of all the components of the bearing. 
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If the force was too large, the cage would accelerate or decel- 
erate until an equilibrium could be reached between all the rollers 
and the cage. 

The model which was used in' the computer program is that given 
by Dowson, Markho and Jones [14] • This model describes the lubri- 
cation of lightly . loaded rigid cylinders by an isoviscous lubricant 
in combined rolling, sliding and normal motion, Fig. 7 . The gov- 
erning equation is the Reynolds equation with side- leakage neglected. 
The boundary conditions used for this model are: 
at the' inlet, the pressure is zero; 
at the outlet, the pressure and the pressure 
gradient are zero. 

The normal pressure distribution is computed and from this the 
normal force and traction forces are determined.* Dowson, Markho 
and Jones presented their results and gave regression equations for 
the instantaneous load carrying capacity and the surface tractions. 
The data given are ideal for inclusion in the roller bearing dynamics 
program because of the relative ease in calculating the regression 
equations • 

The expression for the instantaneous load carrying capacity is 

N = r 1.0 - 1.89351 q + 1.54192 q^ - 0.33529 q^ 

“o 

- 0.17678 q - 2.85050 q H + 1.73361 q^ H 1 (36) 

o oj 

if q ^ 0 
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N == 


4.896 

H 


j^l.O - 1.87439 q + 1.70716 + 0.62039 

+ 0.04271 q - 3,31298 q + 1.35849 q^ hJ ( 37 ) 


if q < 0 


where N is the dimensionless normal load carrying capacity per 

' ' ■ n 


unit length of cylinder = 


liu 


is the dimensionless (normal/entraining) velocity pa- 

.. _ w , ^eq ^l/2 

rameter = ) 


H 


P- 

u 

w 

R 

( 

Q 


is the dimensionless minimum film thickness = 


R 


eq 

n is the normal load per unit length of the cylinder 


is the dynamic viscosity of the lubricant 

U + R n 


is the entraining velocity = 


eq 


eq 


is normal velocity of the cylinder 

is the radius of the geometrically equivalent 

cylinder near the plane, 

is the angular velocity of the cylinder 


The surface traction acting on the roller is given by 

^roller = TW [b ‘ C (38) 




The surface traction acting on the cage pocket surface is 
given by 


cage 


H 




( 39 ) 


29 


where F and F are dimensionless tangential surface forces 

roller cage 


arising from viscous stresses - 


p,u 


B, C 

V 


are constants which are defined below 
is the dimensionless sliding ratio = 

u - R n 


2 (- 




U + R 0 

eq 


■) 


The constants B and C are defined as follows; 

2 


B = [l - 0.76128 q - 0.00358 q^ + 0.42000 q^ 

- 0.21924 q^ - 60.53144 q + 19.56741 q^ 


(40) 


if q ^ 0 

B = [l - 0.76443 q - 0.02328 + 0.37387 q^ 

+ 0.17758 q^ - 60.92918 q - 21.20263 q^ ] (41) 

if q < 0 


where 


(B) „ = 4.5685 - 681.56 H + 341795 H ^ 

' q— 0 o o 

C = (C)q^Q [ 1 - 0.19300 q + 0.08479 q^ + 0.33245 q^ 

- 0.11347 q^ - 4.03120 q H + 5.55690 q^ H 1 
^ o oj 

if q ^ 0 


(42) 


(43) 
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C = (C)q:,0 ■ 0*19648 q + 0.07182 + 0.19631 q^ 

+ 0.08557 q^ - 3.71562 q - 1.15432 q^ (44) 

if q < 0 

where 

(C)q =0 " 3.4843 - 113.66 + 56903 (45) 

These expressions are generally valid for the range of values of the 
variables (q) and (H^) of 


-1.0 < q < + 0.75 


10 


-6 


< H < 10 


-3 


When the variable goes below 10 , the model is switched 

over to a Hertz contact model given by Loo [ 15 ] and linearly approx- 
imated by Conry and Seireg [16] for two cylinders in contact; 


P 



2 (l-v^) 


)(■ 


-J— '^W 

12.5 TH 


(46) 


where P ' 


is the normal force per unit width of the roller 

resulting from the indentation w 

• H 


w is the indentation of the cage by the cylinder 
H 

E is the modulus of elasticity of the roller and the 
cage 

V is Poisson's ratio 
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The indentation of the cage by the cylinder is measured from 
the point where = 1.0 X 10 ^ and the resulting normal force per 
unit width (P' ) is added to the normal force per unit width (n) 
evaluated at = 1.0 X 10 in equations (36) and (37). 

The minimum film thickness and the indentation (if any) of 
the roller relative to the cage in the cage pocket is defined as 


h 

o 



- d 


(47) 


and 


h 



eq 


> 10 





h -l.OxlO^j^ 
o ^ 


eql 


(48) 


h , 

if < 10' 

eq 


where 



is the circumferential clearance in the cage pocket 

is the equivalent radius of the roller and cage pocket 
surface 

is the displacement of the center of the roller 
from the center of the cage pocket. 


The surface tractions acting on the surfaces of the cage pocket 

and the roller are taken to be the total normal force times a Coulomb 

— 6 

friction constant of 0,05 when H < 10 

o 

These relationships are programmed in subroutine DIFFUN [5] with 
equations (36) to (45) programmed in subroutines RLCGF and FRICTF [5]. 
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Drivins force between cage and land 

The roller bearing type considered in this report has a cage 
which rides on two inner race lands. Under full hydrodynamic lubri- 
cation, the relative motion between the cage and the inner race 
land gives rise to a pressure distribution between the inner sur- 
faces of the cage and the land surfaces upon which it rides . The 
expression for the pressure distribution in this oil film is given 
by Kirk and Gunter [ 17] . 

The cage riding on an oil film between it and the inner race 
land can be modeled as an inside out journal bearing. The "journal" 
motion is prescribed while the "bearing" (in this case, the cage) 
moves in a path subject to the Newtonian equations of motion for 
the cage. 

The expression for the forces acting on the cage is based on 
equation 3.42 of fl7], but modified to consider prescribed displace- 
ments and velocities of the inner race. 
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where 


F and F 
X y 

li. 

^ shaft 
u>c 

X and y 
and Ay 
C 


are the hydrodynamic forces on the cage 

is the dynamic viscosity of the lubricant 

is the half -width of the total cage/land surface 

is the angular speed of the inner race 

is the angular speed of the cage 

are the displacements of the cage relative to 

a fixed frame 

are the displacements of the inner race relative 
to a fixed frame 

is the radial clearance between the cage and the 
inner race land 


The integral in equation (49) is evaluated by a composite 
Legendre - Gauss quadrature using a 96 - point formula. This eval- 
uation is performed subject to the restriction that the argument of 
the Integral is negative. If it is positive, it is set to zero over 
that interval for purposes of the numerical integration. This cor- 
responds to integrating over the positive pressure region only. 

These expressions are programmed in subroutines DIFFUN, 
CGPRES and GQU3Z [5], 


Drag torque on the cylindrical roller surface 

As the roller spins about its own axis a shear stress is 
generated on the surface as a result of the shearing of the lubri- 
cant /air mixture present in the open spaces of the bearing. When 
summed over the cylindrical area and multiplied by the value of 
the radius of the roller, the drag torque is obtained. This model 
for the drag torque on the roller is given in Rumbarger et. al. [3] 
and is based on a vortex- turbulent correlation. 
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Before any quantities can be computed, the average clearance 
between the roller and the cage and between the roller and the race 
lands must be determined. Fig. 8. For this computer program the 
clearance between the roller and the cage was left out of this aver- 
aging process as the drag between the roller and cage was considered 
in a previous section. This clearance is defined by equation (50) 
and is the average clearance over the regions I, II, IV and V shown 
in Fig. 8, 


1 

S ’ (2TT - 402) 


r 

j C(0)d© 
circumference 
not adjacent 
to the cage 


The 

average clearance is 

then computed as: 


/ 1 + sin 
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(51) 


where 0^ and 0^ 
and where C 

g 


are defined in Figure 8 

is the total circumferential clearance between 
the roller and the cage pocket as shown in Fig. 


8 . 


The vortex- turbulent correlation was utilized by inserting the 
appropriate values of clearance, viscosity, speed and roller geometry 
into equations (29), (30) w.nd (31) which give the Reynolds number, 
the laminar friction factor and the Taylor number respectively. The 
vortex turbulent correlation is given as [3] . 

0.539474 


= 1.3 




(52) 
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and the effective density of the oil/air mixture about the roller for 
this correlation is 

Peff = 0 (DECFUL) (53) 

If the Reynolds number is less than 2500 or the Taylor number 
is less than 41, the appropriate friction factor is taken to be the 
Iciminar friction factor f (equation (30)). 

Li 

The average shear stress acting on the surface of the roller is 



where Q is the angular velocity of the roller about its own axis. 

The average drag torque on the roller is calculated by integrating 
a differential drag torque over the area of the roller, excluding 
the area between the roller and the cage,. 

W = '■» "" ^ - 402) (55) 

The expressions given in this section are computed in subroutine 
R1245 [5], 

Drag Torque on the roller ends 

Drag torques resist the motion of the roller as it moves through 
the oil/air mixture as was seen in the previous section. In addition 
to drag torques on the cylindrical surface, there are drag torques on 
the roller ends. Following Rumbarger et. al, [3], the effective 
density of the oil/air mixture is given by equation (53). 
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The drag torque is computed by using equation (34) , supported 
by equations (32), (33) and (35). The angular velocity to be used 
is the angular velocity of the roller about its own axis. 

These relationships are programmed in subroutine ROLEND [5]. 

Forces between roller end and guiding flange 

Roller bearings with the roller guided on the inner race subject 
the roller to both driving forces and drag forces due to the close 
clearances and resulting lubrication forces between the roller ends 
and the guiding flange. The flow field in this location is very 
complex but a simplified model is used to approximate these forces 
[3]. 

The forces in this program are based on the shearing of the lub- 
ricant between the roller end and the guiding flange. The roller is 
assumed to move with equal axial clearance on both sides. 

Referring to Figure 9, the area of the roller which projects on 
the guiding flange is divided into 5 horizontal strips and 8 vertical 
strips, symmetrical about the vertical line through the center. 

The effective density of the oil/air mixture is again given by 
equation (53). 

The forces and torques acting on the strips are calculated from 
the shear stress acting on the strip in the direction of the strip 

= Peff^ 
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where -r — is the velocity gradient of the velocity in the direction 

oy 

of the strip. 

The driving and retarding forces and torques are calculated in 
subroutine LANDLIP [5], 

/ 
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THE TWO DIMENSIONAL GOVERNING DIFFERENTIAL 
EQUATIONS OF MOTION 

Consider a roller bearing with (NZ) rollers, guided by a flange 
attached to the inner race and separated by a cage. The outer race 
is assumed to be fixed and the inner race is rotating at some angular 
speed 

For the purposes of writing the two dimensional governing dif- 
ferential equations of motion, the following assumptions are made; 

(1) The motion of the roller is only in the transverse plane 
and is only rotational, both about the axis of the roller and about 
the axis of the bearing. The roller is assumed to be fixed in posi- 
tion between the inner and outer race, and for the case of high 
speed bearings, the roller is always pressed against the outer race 
due to the action of the centrifugal forces. 

(2) The motion of the cage is only in the transverse plane 
and is rotational about its own axis and translational in the plane. 

From these two sets of assumptions, the number of degrees of 
freedom for the roller bearing system is; 

degrees of freedom = 3 + 2(NZ) (57) 

For the purposes of this program, the quantities of interest 
as cage state variables are the angular velocity, angular displace- 
ment, the translational velocities and the translational displace- 
ments. The quantities of interest as roller state variables are 
the angular velocities of the rollers about the roller centers 
and the angular velocities and angular displacements of 
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the center of the roller with respect to the bearing center. This 
gives 3(NZ) + 6 state variables which are required to describe the 
two dimensional motion of the roller bearing. This in turn implies 
that the solution of 3(NZ) + 6 first order differential equations 
are required to provide the state variables which describe the 
motion of the bearing. 

The equations of motion are based on the dynamics of rigid 
bodies. The forces acting on the rigid bodies, the rollers and 
the cage, may be due to contact or deflection of the surfaces 
but the deformations are so small, less than 0,1 percent, that 
the inertial properties are considered to be constant. 

The differential equation of motion of the i*’^ roller 
about its own axis is 

I Cl* Yj 


Referring to Fig, 3, equation (58) becomes 
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i = 1, NZ 


42 



are traction forces calculated by equations (24) 


where 


F and F 
^IN ^OUT 


^TBR ^TFR drag forces on the roller exerted by the 

cage and are calculated by equations (38) or 
a coulomb friction as determined by the sec- 
tion on roller-cage interaction forces, 

Triand ^ driving torque exerted by the guide 

flange on the roller. It is computed in sub- 
routine lANDLiP [5], 


T 


rend 


T 


rcyl 


Ci 


is a drag torque acting on the ends of the 
roller. It is calculated in subroutine 
ROLEND [5]. 

is a drag torque acting on the surface of the 
roller. It is calculated in subroutine R1245 [5] . 

til 

is the angular acceleration of the i 
roller about its own axis 


m 


is the mass of the roller 


r 


is the roller radius 


The differential equation of motion of the i*"^ roller about 
the center of the bearing is, considering the roller as a rigid body 
with the center of mass at the center of the roller: 


m (r + 


R) 


<f>. = V 
/ 1 / 


Torques 


( 60 ) 


where R is the radius of the inner race 
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Referring to Fig. 3, equation (60) becomes 


(r + R)^ = (R + r)^ 


^IN ^OUT 


+ F 1 

B F HLANDJ 


(61) 


1 = 1 , 


NZ 


The differential . equations of motion of the cage describe the 
planar motion of the cage. As the motion of the cage is limited 
to three degrees of freedom, the governing equations, relative to 
a fixed coordinate system, are: 


m X = y F 
cage L X 


(62) 


m 

cage 




(63) 


"cage 




= V 


Torques 


(64) 


where m and I are the mass and polar moment inertia of 

cage cage 

the cage respectively. 

The quantities ^F^ and ^F^ are the sums of forces acting on 
the cage in the x and y directions respectively. The forces 
included in these quantities are the hydrodynamic forces between 
the cage and the inner race land and the forces from the roller- 
cage interaction. The gravity vector is in the positive y direction. 
The forces between the cage and the inner race land are given 
in equation (49). The roller-cage interaction forces acting on the 
cage are calculated from equations (36) through (48) . Referring to 
Figs. (3) and (10); 
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Figure 10. 


The Normal and Frictional Roller-Cage Forces Acting 
on the Cage at the j Pocket 
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( 66 ) 


The torques which act on the cage are classified into three 
categories; (1) the torques due to cage-roller interactions, (2) the 
driving torque due to inner race land-cage interaction, and (3) the 
drag torques about the sides and outer surface of the cage due to 
the shearing of the oil-air mixture. The right hand side of equation 
(64) is written as 
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^ Torques = JJ T dA 

cage- land 
area 

NZ 

+ (R + r) ^ (67) 

j-1 

- T „ - T 

COS CDO 

where T and T are given by equations (25) and (34) respectively. 

I J kj Vf 


The shear stress used in the integral in equation (67) is 
represented as 


_ ^ ^'-“shaft " 

j - 

h 


( 68 ) 


where r^„ 
IN 


h 

M' 

^shaft 




is the inner radias of the cage which rides on 
the land 

is the film thickness 

is the absolute viscosity of the oil 

is the angular speed of the shaft 

is the angular speed of the cage 


Letting 


h = C (1 + e cos ©) 

the integral in equation (67) becomes: 
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where C is the radial clearance between the cage and the land 

L is the half width of the total cage- land surface 

g is the ratio of the minimum film thickness to the 
radial clearance 

Some of the differential equations described above are reform- 
ulated for a more efficient numerical solution. The motion of the 
center of mass of a roller is necessarily limited by the constraints 
of the cage pocket. A critical quantity used in the determination 
of the cage-roller interaction forces is the relative distance be- 
tween a roller surface and its cage pocket surfaces. This quantity 
is obtained by subtracting the angular displacement of the mass 
center of the roller from the angular displacement of the cage, 
which numbers can become very large and which can lead to errors 
in accuracy. 

Rather than solve equations (61) and (64) and subtract the^^'s 
from/i^, equations (61) are replaced by the difference between the 
angular accelerations of the rollers and the angular acceleration 
of the cage. When integrated twice this yields an angular displace- 
ment difference. The calculation will be stable, the displacement 
will be bounded and the numerical problem of differencing large 
numbers will be avoided. The new set of differential equations ob- 
tained from equations (61), (64), (67) and (70) are 
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cage cage 


i = 1, . . . , NZ 

The governing differential equations for the planar motion of 
the roller bearing system are now given by equations (59) , (62) , 
(63), (64) and (71). These are 2(NZ) + 3 equations which corre- 
spond to the number of degrees of freedom given in equation (57), 

Before attempting a numerical solution of this system of equa- 
tions, they will be made non-dimensional. The forces will be 
scaled by a quantity F* which is taken to be the maximum force on 
any of the rollers [2]. This maximum force is related to the mass 
of the roller, a characteristic dimension, in this case, R +r, 
and a characteristic value of time by 

m (R + r) ( 1 )^^ = F* (72) 

where is a characteristic angular velocity 

The dimensionless time T is related to real time t by 


49 


(73) 


T 




Equation (71) is related to the relative displacement of the 
center of the roller in the cage pocket by multiplying it by (R + r) 
and scaling it with the total clearance between the roller and the 
cage pocket (Cg). The relative dimensionless displacement is de- 
fined as 


(R + r) 




Ai = 


(74) 


In addition to scaling forces by F*, torques are scaled by F* r 
and the dimensionless torques and forces are denoted by a supra 
bar (-). 

Equation (59) becomes 
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Equation (71) becomes 
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Define the dimensionless displacements of the cage as: 
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Equations (62) and (63) are then given as; 
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The differential equation relating the angular acceleration of 


the cage to the torques acting on the cage (equation (64)) is 
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The system of differential equations (equations (75), (76), (78) 
(79) and (80))are then solved numerically to obtain the state 
variables of interest. Equations (75) are only integrated once 
to find the angular velocity of the roller. The angular dis- 
placement of the roller about its own center is not needed for 
this problem. Equations (76) are integrated twice to find the 
relative velocity and relative displacement of the roller in the 
pocket. Equations (78), (79) and (80) are also integrated twice 
to find the velocities and displacements of the cage. 
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THE THREE DIMENSIONAL GOVERNING DIFFERENTIAL 
EQUATIONS OF MOTION 

The equations of motion to be presented in this chapter build 
on the governing equations developed in THE TWO DIMENSIONAL GOVERNING 
DIFFERENTIAL EQUATIONS OF MOTION. Four more degrees of freedom are 
added to each roller, an axial translation, a radial translation, a 
skewing of the roller relative to the inner race and a tilt of the 
roller relative to the inner race. 

The cage motion is assumed to remain planar. The cage in a 
typical very high speed bearing has a large angular momentum vector 
which would require torques of high magnitude to force the cage into 
a significant precession. Kingsbury [18] observed this phenomenon 
of planar cage motion in tests of bearings for gyroscopes. 

The planar model of the cage is also used because of the lack 
of good models to relate the potentially skewed surfaces between 
roller and cage to torques about the transverse axes of the cage. 

The roller is so constrained by the side flanges that, for the pur- 
pose of calculating the squeeze film forces, it remains parallel to 
the cage pocket surface. 

In the previous section, the roller motion in the radial di- 
rection was constrained to satisfy the static equations of equi- 
librium of the roller bearing system. In this section, that con- 
straint is relaxed to allow accelerations of the roller in the 
radial direction because of the cross couplings of the tilt angle 
of the roller, the radial deflection of the roller and the normal 
force and tilt torque applied to the roller. These quantities are 
needed to solve the coupled Euler equations of motion for each 
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roller. The static equations of equilibrium are used to establish 
the initial conditions for the roller radial displacements of the 
centers of mass of the rollers and the roll angles of the rollers. 

The complete motion of each roller is obtained by considering 
the motion of the center of mass of the roller in the inertial 
coordinate frame (p , ^ , z) and the rotational motion of the roller 
about its center of mass in the roller frame (^, TIj under the 

assumption that the tilt angle ^ and the skew angle T] are very 
small. These frames are selected such that the p vector and the T] 
axis are parallel and the ^ axis is parallel to the inertial axis z. 

Denoting the mass of the roller by m, the polar moment of 
inertia of the roller (about the ^ axis) by I and the transverse 
moment of inertia of the roller by the equations of motion of 
each roller are [19] ; 
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is the orbital velocity of the roller 


where 




G , G and G as the torques applied about the (F, 'n» 0 axis 
I r\ ^ ^ 

The third of equations (82) are identical to equations (58) . The 

torques G^^^ are given in the right side of equation (59). 

The axial motion of the roller relative to the inner race is deter- 
mined using a kinematic relationship. The axial velocity of the center 
of mass of the roller relative to the race is equal to the product of 
the relative velocity of the center of mass of the roller and the race 
in the Tj direction (U^^^) and the angle of skew of the roller (T]^)» thus 


dz . 
1 

dt 


■ ‘'rel ’ll 


(83) 


i = 1 NZ 


The 


radius vector in equations (81) is given as 


= r + R + 6^^^ 
r 


(84) 


where 


^ is the deflection of the roller mass center and r and 
r 

R are the roller radius and inner race radius respect- 
ively. 


The forces are determined by substituting the 5 ^^^ for 5 J. 

p ^ 

and the Tlj^ for 3 j fo equations (4) , (5) , and (10) with = 0. The 
residual F . is interpreted as the net applied force on the roller 
in the radial direction. The torques, are determined by evalu- 
ating equation ( 11 ) where the residual, F^ + j;j 2 » negative of 

the net tilt torque on the roller. 
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The dimensionless form of the forces, are given in the 

brackets on the right hand side of equations (61). If were a 
constant, the second of equations (81) would be identical to 
equations (61). 

The torques in equations (82) have two components, a 

torque due to unevenly distributed traction forces along the length 
of the roller at the roller-race interfaces and a torque due to 
elastic deformation between the edges of the roller and the guiding 
flanges of the inner race when the roller impacts those guide flanges. 
The torque due to unevenly distributed traction forces is continuous 
in time. The misalignment of the bearing inner race to the outer 
race contributes to the uneven load distribution on the roller which 
in turn affects the magnitude of the traction coefficient. The net 
moment of the product of the normal loads on the "slices" of the 
roller (from Chapter 2) and the traction coefficients for the same 
slices is the torque which tends to skew the roller. 

The torque due to the impact of the roller edge(s) with the 
guide flange (s) gives rise to a restoring torque which can have a 
large magnitude. The basis for calculating the contact force is 
given by Brewe and Hamrock [20] where an approximate force relation- 
ship is given as a function of the geometry of the contacting sur- 
faces and the material constants. 
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where f = 1.5277 + 0.6023 log (R /R ) 
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where 



Ax 


k 


f 

E’ 


is the contact force due to impact 
is the contact deflection 

is an approximation of an elliptic integral of 
the second kind 

is the effective radius 
etc. are radii of curvature 
is the ellipticity parameter 
is Poisson's ratio 

is an approximation of an elliptic integral of 
the first kind 

is the effective modulus of elasticity 


The subscripts A and B refer to the two bodies in contact, respectively. 
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The equivalent radius in one principal plane at the contact 
point is approximated as; 



(94) 


The equivalent radius in the other principal plane at the 
contact point is approximated as: 



(95) 


Therefore 


from equations (89), 

JL ^ _1 

^eff ^fm 


(94) and (95) 
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where r is the minor radius at the corner of the roller 

m 

r is the minor radius of the inner race guide flange 

fm 

at the corner. 


When the roller is skewed, the contact between the roller and 
the guide flange is at a point located at the intersection of the 
roller circumference and the circumference of outer surface of the 
guide flange. 

The restoring torque after impacting the guide flange is 

Tf = Fj (97) 


as shown in figures (4) and (11). Figure (11) also illustrates that 

(98) 


r_ = r sin 
f F 


0„ = arc cos 
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2 2 2 
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2r (R + r) 


where R^ 


In summary 


is the radius of the outer surface of the inner race 
guide flange 

, the torques may be expressed as; 
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Figure 11. The Geometry of the Roller Bearing System Used to 
Establish the Point of Contact Between a Skewed 
Roller and the Guide Flange 
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where 
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,(ij) 
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and F 


(ij) 

OUT 



is the distance from the center of mass 
of the roller to the center of the "slice" 
as described in THE CONSTITUENT MODELS, 
are the tractive forces calculated over 
the contact area of the j ‘th slice using 
equations (24) 

is the contact force between a roller 
corner and the side walls of the guide 
flange. 


In equation (99) the four potential contact points of the roller 
are considered. Fig. 4. In the program, the locations of these 
points are monitored and if contact occurs, the corresponding 
is calculated. If there is no contact, the corresponding is set 
to zero. 

From the discussion presented thus far in this chapter, ten 
first order differential equations are required to describe the motion 
of each roller. For the bearing with NZ rollers and a cage in planar 
motion only, the total number of first order differential equations is 
6+10 (NZ). 

The remainder of this section deals with the simplifications 
made to this system of equations. The essential character of the 
motion is retained while ignoring some smaller effects. 
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Consider now an harmonic solution of the first of equations 


(81), after the substitution of equation (84) to be 

# 

5^^^ = cos t 


( 100 ) 


Then 

m ( - Di cos t - (R + r) 4 ^^ cos t) = (101) 


The quantity on the left hand side is dominated by the term (R + r) , 

since 


R + r » (102) 

-3 -2 

The value of would typically be of the order of 10 or 10 and 

R + r is of the order 1. This implies that the forces should be 

P 

equal to the centrifugal forces, the condition enforced in the static 
load distribution problem considered in THE CONSTITUENT MODELS. Thus 
the first of equations (81) may be eliminated. 

At each point in time every roller is placed in static equilib- 
rium which implies that the velocity and acceleration of the center 
of mass of the roller are zero in the radial direction. Thus, the 
second of equations (81) become, using equation (102) 


m (r + R) 



(103) 


i = 1, . . . , NZ 


Equations (103) are identical to equations (61). 
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The first of equations (82) are eliminated because of the enforce- 
ment of the static equations of equilibrium. The torques are zero 
as a result of this constraint and the values of the roll velocity and 
acceleration of the rollers are taken to be zero. The second term 
of the first of equations (82), T]^ is estimated to have a mag- 

_ 9 _ 2 

nitude of the order of 1.13 X 10 N-m (10 in - ^b^) which would entail 
only a minute change in the static normal load distribution of the rollers. 

The second term on the left hand side of the second of equations 
(82) is set to zero because of assumption of zero tilt velocity and 
acceleration of the rollers. This set of differential equations for 
the skew displacements becomes: 


Tii 


T| 


(104) 


i = 1, . . . , NZ 


The simplified system of governing differential equations for 
the three dimensional motion of the rollers is given by the third of 
equations (82), equations (83), (103) and (104). The equations of 
motion of the cage remain unchanged from the planar motion case. 

The dimensionless form of equations (83) is then written as 
using equations (73) and (74) 


d z 


dT 


= -[ 






(R + r) 


R 


^^shaft 

C 

z 



(105) 


where z . - -~ 
1 C_ 


62 



is the total axial clearance between the roller ends and 


C 

z 


^shaft 


the inner race guide flanges 
is the shaft angular velocity 


The dimensionless form of equations (104) for the skew displace 
ments are of the form 


dT'^ 


F* 


O) , 


NP 

Y 
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^(j) 

^ IN OUT ^ 


+ 


F* 




(10i3) 


The reduced system of dimensionless differential equations to 
be solved is equations (75), (76), (78), (79), (80), (105) and (106) 
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DISCUSSION OF RESULTS 


The equations of motion developed in this report differ in certain 
respects from those presented by Gupta [4] . Gupta presents the three 
equations of motion of the roller mass center in an inertial reference 
frame and presents the three Euler equations of motion in a roller fixed 
frame. He then presents a transformation matrix to convert the moments 
from an inertial frame to a roller fixed frame. His formulation is gen- 
eralized and contains the full six degrees of freedom for each roller. 

Two formulations of the equations of motion were presented in this 
report, the two dimensional problem and the three dimensional problem. 

The motion of the cage is limited to two dimensions for both cases. The 
two dimensional formulation only permits planar motion of the roller. The 
radial motion is neglected as the races and centrifugal forces effectively 
constrain it to be very small. 

The three dimensional formulation again neglects the radial displace- 
ment and the tilt angle dynamics of the roller. The tilt of the roller 
is assumed to be a function of the misalignment of the races, the loading, 
and the geometry and elasticity of the rollers. The axial motion of the 
roller relative to the race is expressed by a kinematic condition rather 
than a dynamic condition. 

The equations of motion which describe the motion of the roller center 
in the plane of the bearing are modified in this report. Many of the forces 
which act on the cage and the rollers are functions of the relative displace- 
ments between roller and cage. To keep numerical stability, new differential 
equations are formulated in terms of these relative displacements. Gupta 
[4] makes no mention of using this approach. 
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The approach used by Gupta [4] for calculating roller-race inter- 
action is essentially the same as used in this report. The work of 
Palragren [8] was used as the basis for the normal load distribution cal- 
culation. The lubricant model used in this report was for Mil-L-7808 
type lubricant while Gupta [4] used Mil-L-7808 and 5P4E polyphenyl ether. 

The roller end to race flange interactions are treated in similar 
fashion in this report and in Ref. [4]. The roller to cage and race to 
cage interactions are also treated in a similar fashion to Ref. [4]. The 
churning and drag loss models in this report and in [4] are based on the 
work of Rumbarger [3] . These are approximate models which depend on an 
effective density. This approach has been used in most rolling bearing 
applications. The effective oil density is in effect a parameter which 
must be adjusted to the operating conditions which are being simulated. 
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CONCLUDING REMARKS 


The two and three dimensional differential equations of motion of 
a high speed cylindrical roller bearing have been derived in this re- 
port. Simplifying assumptions were made which were specific to the 
cylindrical roller bearing. The original differential equations were 
manipulated to provide a more numerically stable system of differential 
equations . 

Approximate models were utilized to compute churning losses, traction 
forces, roller to cage and roller to race interaction forces. As more 
knowledge becomes available in these specific areas, changes to the over- 
all differential equations can be made. 
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h 

h 

o 

k 

Z 

m 

m 

cage 

n 
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radial clearance (m) 

distance between edge of cage pocket to center of 
roller (m) 

approximation of an elliptic integral of the second 
kind 

traction force (N) 
friction factor 

2 

acceleration of gravity (m/s ) 

total contact deformation at the inner contact of the 

i^^ slice of the roller at angular location (m) 

total contact deformation at the outer contact of the 

i^^ slice of the roller at angular location ^ (m) 

film thickness (m) 

minimum film thickness (m) 

ellipticity parameter 

length of a roller (m) 

mass of a roller (kg) 

mass of the cage (kg) 

normal load per unit length (N/m) 

dimensionless velocity parameter 

load intensity (N/m) 

2 

maximum Hertzian contact pressure (N/ra ) 
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u 
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V 
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H 
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roller radius (m) 
characteristic radius (m) 

distance as described in equation (98) (m) 
minor radius of the inner race guide flange at the 
corner (m) 

inner radius of cage (m) 

minor radius at the comer of the roller (m) 
outer radius of cage (ra) 

mass average velocity of fluid, equation 26 (m/s) 
entraining velocity (m/s) 

sliding velocity at the roller- race contact (m/s) 

normalizing constant for sliding velocity (m/s) 

rolling velocity at the roller-race contact (m/s) 

dimensionless sliding ratio 

width of a slice of the roller (m) 

normal velocity (m/s) 

normal indentation of a roller (m) 

absolute displacement of cage in x-direction (m) 

absolute displacement of cage in y-direction (m) 

location of the i slice of the roller (m) 
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F 

cage 


F 


F 


^HLAND 
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constant equal to 5.46 X 10 , equation (1) 

surface area (m^) 
constant, equations (40) and (41) 
constant, equations (42) and (43) 
radial clearance (m) 

circumferential clearance in the cage pocket (m) 
average clearance between roller and cage (m) 
initial separation due to crowning of the slice 

of a roller (m) 

friction factor, equation (35) 

axial clearance between roller ends and the inner 
race guide flanges (m) 
an arbitrary constant 

ratio of the oil volume in the bearing to the total 
volume of the bearing 

2 

modulus of elasticity (N/m ) 

2 

effective modulus of elasticity (N/m ) 
maximum force on the rollers, used for scaling (N) 

normal forces on the roller exerted by the cage (N) 

centrifugal force on a roller (N) 

dimensionless traction force acting on cage pocket 
surface 

normal forces on the roller exerted by the cage (N) 

force exerted on the roller sides by the inner race 

land (N) 
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contact force due to impact (N) 

traction force between roller and inner race (N) 
residual quantity used in equations (9) and (10) 
traction force between roller and outer race (N) 
dimensionless traction force acting on roller 
drag forces on the roller exerted by the cage (N) 
drag forces on the roller exerted by the cage (N) 
hydrodynamic force on cage in x-direction (N) 
hydrodynamic force on cage in y-direction (N) 
torque (N-m) 

dimensionless minimum film thickness 

mass moment of inertia of roller about its axis of 
2 

rotation (kg - m ) 

2 

mass moment of inertia of cage (kg - m ) 

2 

transverse moment of inertia of roller (kg - m ) 

half-width of total cage/land surface (m) 

dimensionless normal load carrying capacity 

til 

total normal force on the j roller at the inner 
race contact (N) 

til 

total normal force on the j roller at the outer 
race contact (N) 

number of slices along the axis of roller 
Reynolds number 
Taylor number 

number of rollers in bearing 
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R 
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T 
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^CDS 

^rcyl 

"^rend 

T 
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normal force per unit width (N/m) 
diametral clearance of the bearing (ra) 
applied load to a roller (N) 
radius of the inner race (m) 
effective radius 
equivalent radius (m) 

radius of the outer surface of the inner race guide 
flange (m) 
dimensionless time 

drag torque on outer cylindrical surface of cage (N - m) 
drag torque on the sides of the cage (N - m) 
restoring torque (N - m) 

drag torque acting on the surface of the roller (N •- m) 
drag torque acting on the ends of the roller (N -m) 
driving torque exerted by the guide flange on the 
roller (N - m) 
surface velocity (m/s) 
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T1 

T1 

P 
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V 

V 

I 
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P 

P 

^eff 

T 
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one-half of the included angle of the cage pocket 
(radian) 

tilt angle of the j roller (radian) 
contact deformation (m) 

relative approach between the roller axis and the 
outer ring along the radial line passing through 
the center of the roller at angular location 

j 

eccentricity ratio 
principal axis in roller frame 
angle of skew (radian) 

angular displacement of the cage (radian) 
dynamic viscosity of oil (N-s/m ) 
traction coefficient 

normalizing constant for traction coefficient 
2 

oil viscosity (m /s) 

Poisson's ratio 

principal axis in roller frame 
angle of tilt (radian) 

radial direction of roller motion in an inertial frame 
displacement of roller in the radial 
direction (m) 

3 

density (kg/m ) 

3 

effective density of oil/air mixture (kg/m ) 

2 

shear stress (N/m ) 

2 

wall shear stress (N/m ) 
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^shaft 

r 

Ax 


e 

X 
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angular position of the j roller relative to the 
x-axis (radian) 

angular velocity of the i^^ roller about its own 
axis (m/s) 

principal axis in the roller frame 
angular velocity 

angular velocity of the cage (radian/s) 

characteristic angular velocity, equation (72) (radian/s) 
angular velocity of the shaft (radian/s) 
constant, equation (11) 

component of radial displacement of bearing inner race 
in x-direction (m) 

component of radial displacement of bearing inner race 
in y-direction (m) 

angular misalignment of bearing inner race relative to 

the outer race about the x-axis (radian) 

angular misalignment of bearing inner race relative to 

the outer race about the y-axis (radian) 

angular location of an axis with respect to the x-axis 

about which the total angular misalignment is applied 

(radian) 

approximation of an elliptic integral of the first kind 
angular velocity (radian/s) 
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Subscripts 


L 

C 

C 

n 

p 

< 1 > 


laminar 
5 direction 
^ direction 
fl direction 
p direction 
(j) direction 


z 


z direction 
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